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Abstract 

We use latest data from solar system planetary orbital motions to put constraints on 
some Galileon-induced prcccssional effects. Due to the Vainshtcin mechanism, the Galilcon- 
type spherically symmetric field of a monopolc induces a small, screened correction oc y/r 
to its usual r^ 1 Newtonian potential which causes a secular precession of the pericenter of 
a test particle. In the case of our solar system, latest data from Mars allow to constrain 
the magnitude of such an interaction down to a < 0.3 level, where a corresponds to the 
non minimal coupling of the Galileon to matter. Another Galileon-type effect which might 
impact solar system dynamics is due to an unscreened constant gradient induced by the 
peculiar motion of the Galaxy. The magnitude of such an effect, depending on the different 
gravitational binding energies of the Sun and the planets, taken into account by £, is £ < 
0.004 from the latest bounds on the supplementary perihelion precession of Saturn. 

PACS: 04.80.-y; 04.80.Cc; 04.50.Kd; 95.10.Km 

1 Introduction 

Galileon theorjQ [2] , originally formulated in 4D flat spacetime and later generalized also 
to curved backgrounds [3], yields a modified gravity model relying upon arjl additional light 
scalar field ir, the Galileon, endowed with derivatve self-interactions^|; for a recent overview of 
this— and of many others— modified model of gravity, see, e.g., [5]. The name comes from the 
fact that the Lagrangian for ir is left unchanged by a transformation of it which generalizes the 
usual Galilean invariance. The Galileon scenario aims to explair0 the observed acceleration of 
the Universe [8j El \W\ \TT\ [T2| [13] without resorting to Dark Energy by modifying the known 
laws of gravitational interaction at large distances. At the same time, it avoids to compromise 

1 See, e.g., pQ for a recent review. 

2 At least, in its original formulation. Multi-Galileon scenarios, introduced for the first time in [4], appeared 
to cope with certain potential inconsistencies of the single- Galileon model. 

They can help in preserving the observed agreement of the General Theory of Relativity with solar system 
phenomenology by partly screening the effects of the extra-scalar field at relatively short distances. See below 
and [5]- 

4 See [51(7] for an earlier formulation in a different context. 
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their agreement with observations at solar system scales through a screening based on an 
implementation of the Vainshtein mechanism [J3] . The Vainshtein mechanism was first shown 
to work in all its generality [TS], and then in the decoupling limit as well |16j . Inspired to 
the multidimensional braneworld model by Dvali, Gabadadze and Porrati (DGP) \17\ fT8j [T9"l 
I2(J| . the Galileon theory, which has been shown the emerge explicitly from higher-dimensional 
braneworld scenarios and massive gravity [211 [22l [23] , is able to cure certain drawbacks of the 
latter like the appearance of ghosts in the self-accelerated branch [25]. Nonetheless, the single- 
Galileon scenario is not entirely consistent because of certain issues appearing both at classical 
and quantum scales, along with unwanted superluminal features [5]. A bi-Galileon theory [25] 
seems able to successfully face such problems. 

In this paper, we want to effectively put quantitative constraints on some consequences 
of the Galileon scenario with orbital motions in our solar system. For other constraints on 
Galileon-type interactions inferred from cosmological observations of different phenomena, see, 
e.g., [2B1 [27] [2S]. As an exciting perspective, local observations can fruitfully shed light on 
the Dark Energy problem and nature of gravity. Indeed, as we will show in the next Sections, 
Galileon induces various kinds of local, non-negligible orbital effects, given the present-day 
level of accuracy in constraining non-standard planetary orbital precessions from solar system 
observations |29j. In particular, in Section [2] we look at the small orbital precessions induced 
by a single Galileon-type field of a central monopole. In Section [3] we inspect the effects caused 
by the unscreened Galileon constant gradient due to the large scale-induced peculiar motion 
of the Galaxy. As a result, a small differential Sun-planets extra-acceleration should occur in 
view of their different gravitational self-energies. For other proposed tests involving different 
astronomical systems like compact objects in the Milky Way crossing the Galactic plane, and 
stars and the central black hole in M87, see [30]. As a general cautionary remark in using 
certain scenarios less directly accessible and known with respect to our solar system, we want 
to note that they may typically suffer from large systematic effects whose accurate knowledge 
is often lacking, contrary to Sun's planetary arena. Moreover, in some cases they are also quite 
model-dependent in the sense that they heavily rely upon theoretical assumptions which are 
still speculative since they have not yet been tested independently with a variety of different 
phenomena, or have not yet been directly tested at all. It should also be remarked that some 
astrophysical phenomena like the behavior of black holes in extragalctic scenarios may crucially 
depend on the composition, formation and dynamical history of the systems considered. This 
is why we adopt well known and largely tested orbital motions of some natural bodies in our 
solar system. Section H] summarize our findings. 

2 The orbital precessions due to a single Galileon-type field of 
a central monopole 

As a result of the Vainshtein mechanism, one of the forms in which the Galileon extra-potential 
for a localized source of mass M can appear is [2 [3TJ [5] 

U Ga x = aH VGMr, (1) 

where G is the Newtonian constant of gravitation and [32] Hq = 73.8 km s -1 Mpc -1 = 2.4 x 
10~ 18 s" 1 is the current value of the Hubble parameter. We emphasize that, in principle, a 
free parameter-a mass scale-which does not necessarily coincide with Hq can be present in eq. 
(PQ). Thus, the final results of our analysis should depend on such a free parameter and the 
coupling constant a, which is expected to be of the order of unity if the aforementioned free 
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mass parameter does not exactly coincide with Hq. eq. ([I]) comes from one of the non- linear 
terms in the spherically symmetric, static equation of motion for tt near an object of mass M 
[2"| 1314 15]: such non-linearities become dominant just at relatively short distances. Notably, it 
is the same non-linear interaction which allows for both the Vainshtein screening and for the 
self-acceleration of the Universe [2J. 

As far as the solar system is concerned, the Galileon becomes dominant with respect to the 
usual Newtonian potential at distances of the order of the Vainshtein radius for the Sun 
which amounts to about 240 pc. Instead, throughout the planetary region it is 
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the screening of the Galileon is, thus, quite effective when the potentials are compared. Nonethe- 
less, it turns out that eq. (pQ) affects planetary orbital motions with secular precessional effects 
whose expected magnitude may not be hopelessly negligible, being, thus, possible to effec- 
tively constrain it with present-day observations. Moreover, the expected improvements in the 
knowledge of the orbital motions of some of the planets of the solar system currently orbited by 
accurately tracked spacecrafts like MESSENGER (Mercury) and Cassini (Saturn) will further 
strengthen the bounds on the Galileon-type interaction. 

In view of the screening of eq. ([2]), the Galileon-induced planetary orbital effects can be 
worked out with standard perturbative techniques [33]. 

The average of eq. ([I]) over a full orbital period of a test particle, obtained by using the 
true anomaly / as fast variable of integration, is 
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wher^l E is the complete elliptic integral and K is the complete elliptic integral of the first 
kind; a and e are the semimajor axis and the eccentricity of the test particle's orbit. 

By using eq. ([3]) in the Lagrange planetary equations [33] for the variation of the osculating 
Keplerian orbital elements of a test particle, it follows 
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(6) 



5 Here tv is not the Galileon scalar field, being, instead, the usual ratio of circumference to diameter in Euclidean 
geometry. 
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where J is the inclination of the orbital plane to the reference {x, y\ plane adopted, Q is 
the longitude of the ascending nod^l, w is the longitude of pericenteiQ and M. is the mean 
anomaljo 

We remark that eq. (jl])-eq. ([9]) are exact in the sense that no a-priori simplifying assump- 
tions on the orbital geometry of the test particle were assumed. Notice also that eq. (jl])-eq. 
depend on the orbital configuration of the test particle through the eccentricity e; they are 
independent of the semimajor axis a. To order 0(e 2 ), eq. ([8])-eq. (J9|) reduce to 
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We notice that the spherical symmetry of eq. (JTJ) would immediately allow to infer that no 
secular effects at all can occur for a,e,I,f2. Indeed, the Lagrange rate equation for a |33j 
contains the partial derivative of the averaged perturbing potential (U peT t) with respect to Ai, 
which is proportional to t. The Lagrange equations for e, /, fl |33j are formed with the partial 
derivatives of ([/pert) with respect to Q,uj,I (and M. as well in the Lagrange equation for e), 
which are absent in all spherically symmetric perturbing potentials. 
From latest Mars data, summarized in Table [H it follows 



a 



0.07 + 0.3 



(12) 



6 It is an angle in the reference {x, y} plane counted from the reference x direction to the line of the nodes, 
which is the intersection of the orbital plane to the {x, y} plane itself. 

7 It is a "dogleg" angle since it is denned as zu = Q + oj, where cj is the argument of pericenter lying in the 
orbital plane and counted from the line of the nodes to the point of closest approach. 

8 It is defined as M = rib (t — to) , where rib = V GM/a 3 is the Keplerian mean motion to is the time of 
passage at pericenter. 
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Table 1: Supplementary precessions Am, AQ of perihelia and nodes of some planets of our solar 
system estimated by Fienga et al. [29] with the INPOPlOa ephemerides. Data from Messenger 
and Cassini were used. All standard Newtonian/Einsteinian dynamics was modeled [29], with 
the exception of the solar Lense-Thirring effect. However, it is relevant only for Mercury, given 
the uncertainties released. The reference {x,y} plane is the mean Earth's equator at J2000.0. 
mas cty -1 stands for milliarcseconds per century. 





A ft (mas cty 1 ) 


Am (mas cty 1 ) 


Mercury 


1.4 ±1.8 


0.4 ±0.6 


Venus 


0.2 ± 1.5 


0.2 ± 1.5 


Earth 


0.0 ±0.9 


-0.2 ±0.9 


Mars 


-0.05 ±0.13 


-0.04 ±0.15 


Jupiter 


-40 ± 42 


-41 ± 42 


Saturn 


-0.1 ±0.4 


0.15 ±0.65 



3 Differential orbital precessions induced by unscreened large 
scale structure Galileon 

According to Hui and Nicolis [3D], the Galileon symmetry allows for the existence of an addi- 
tional constant gradient in any solution of the Galileon equation. In the case of a background 
external field, such an extrslfl V(/2 cx t would impart an extra-acceleration to a body of mass M 



A = -a Ij^j V^xt, (13) 

where Q is a conserved scalar charge which coincides with M if the self-gravitational binding 
energy is negligible; for a compact object like a neutron star Q/M < 1, while Q/M = for a 
black hole. If, for a given object, such a constant gradient actually exists or not depends on 
the boundary conditions of the specific scenario considered. In the case of a galaxy, they are 
yielded by the surrounding large scale structure. Hui and Nicolis [30] pointed out that it would 
be able to induce a long-wavelength scalar field which, on the scale of a galaxy, would resemble 
just a constant gradient. It would be able to penetrate the Vainshtein zone of the galaxy, acting 
unsuppressed on it and on all its constituents which, thus, would experience a differential fall 
according to their Q/M ratios. 

In principle, also a Sun-planet pair should experience such a small differential tug because, 
although non-relativistic, our star has a gravitational binding self-energy much larger than the 
planets. Thus, it is possible to analyze perturbatively its effects on the orbital motion in detail. 
It can be anticipated that, in view of the anisotrpic character of the interaction mediated by 
V(/3 C xt> the resulting orbital perturbations should likely not be limited to just the perihelion 
and the mean anomaly. Indeed, the averaged perturbed potential of the relative Sun-planet 
motion due to the Galileon-induced large scale structure effect is 

(Uiss) = (CV99 cx t • r) = --£aeA jcos u (k x cos Q + k y sin Q^j + 



In [30], it is 7r = cup. 
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where £ accounts for the Sun-planet difference in their Q/M, and 

A = | Vc^extl , 



(14) 
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As expected, the Lagrange planetary equations, applied to eq. (|14|) . yield the following non-zero 
secular rates of changes 
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(22) 



Also the long-term rates of change of eq. (|17p -eq. (|22|) are exact in the sense that no a-priori 
simplifying assumptions on the orbital configuration of the test particle were assumed. They 
are proportional to \fa, so that they are larger for the outer planets. For a more realistic 
contact with actual observations, which are always referred to a specific reference frame, we 
also did not a-priori align the unit vector k of the external gradient to any specific direction. It 
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is important to notice that our results eq. (fT7]) -eq. (|22|) are quite general, so that they do not 
necessarily refer to a Sun-planet pair, being valid for other systems including a compact object 
as well. 

The peculiar motion of the Milky Way [33] can be assumed coincident with that of the Local 
Group with respect to the Cosmic Microwave Background (CMB) rest frame. It occurs at 

VhG = 627 km s" 1 (23) 

towards 

RA = 11.11 hr, (24) 



DEC = -27.33 deg, (25) 

corresponding to 

k x = 0.8717, (26) 

k y = 0.1711, (27) 

k z = -0.459. (28) 

The magnitude of the constant gradient can be approximated to the product of the speed of 
the peculiar motion times the Hubble parameter |30j. In the case of the Milky Way, from eq. 
([23]) we have 

A ~ v LG H = 1.5 x 10~ 12 m s~ 2 . (29) 

From eq. (|2ip and the uncertainty for the supplementary perihelion precession of Saturn in 
Table Q] it turns out 

|£| < 0.004. (30) 



4 Summary and conclusions 

As far as the Galileon-type spherically symmetric field of a central monopole is concerned, we 
looked at the term oc Hq^/t arising from the Vainshtein mechanism. For the Sun and its planets, 
it is screened to a high level with respect to the usual r _1 Newtonian potential; nonetheless, it 
induces non-vanishing secular precessions of the longitude of perihelion and the mean anomaly 
which depend only on the planetary eccentricities. Interestingly, their expected magnitudes are 
rather close to the present-day level of accuracy in constraining the planetary supplementary 
precessions, of the order of < 1 milliarcseconds per century. Thus, we were able to effectively 
constrain the strength parameter a of the residual Galileon interaction down to \a\ < 0.3 from 
Mars data. 

Then, we looked at the unscreened Galileon-type local orbital effects caused by the Galactic 
peculiar motion due to large scale structures. They are proportional to the difference of the 
gravitational self-energies between the Sun and its planets. Because of the anisotropy of such 
an interaction, all the osculating Keplerian orbital elements of a test particle undergo non-zero 
long-term variations, apart from the semimajor axis. By using the supplementary perihelion 
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precession of Saturn, we were able to infer |£| < 0.004 for the strength parameter of such an 
interaction. 

Further gathering and processing data from ongoing and forthcoming spacecraft-based mis- 
sions orbiting some planets of the solar system like Mercury (MESSENGER, BepiColombo) and 
Saturn (Cassini) will allow to enhance such thrilling opportunity of shedding light on crucial 
cosmological issues from local observations. 
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